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Strong influence of water molecules on the transport and magnetic properties of DNA, observed in 
this study, opens up real opportunities for application of DNA in molecular electronics. Interaction 
of the nucleobases with water molecules leads to breaking of some of the tt bonds and appearance of 
unbound tt electrons. These unbound electrons contribute significantly to the charge transfer at room 
temperature by up to 10^ times, but at low temperature the efficiency of charge transfer is determined 
by the spin interaction of the two unbound electrons located on the intrastrand nucleobases. The 
charge exchange between the nucleobases is allowed only when the spins of unbound electrons are 
antiparallel. Therefore, the conductance of the DNA molecule can be controlled by a magnetic 
field. That effect has potential applications for developing a nanoscale spintronic device based on 
the DNA molecule, where efficiency of spin interaction will be determined by the DNA sequence. 



I. INTRODUCTION 



The growing interest on the application of or- 
ganic molecules for building nanoelectronic devices is 
motivated by several factors: conductance proper- 
ties, self-assembly and molecular recognition properties. 
Self-assembly of molecular building blocks into well- 
structured systems allows us to exclude physical manip- 
ulation during fabrication of nanosize devices. Moreover, 
the structural damages of the electrical elements can be 
repaired without any physical contact using the molecu- 
lar recognition. Since the DNA molecules satisfy all these 
requirements and, moreover, have been found to conduct 
charge, its application to developing a nanoscale molec- 
ular devices is quite promising. While pristine DNA is 
not a good conductor, its conductivity can be greatly en- 
hanced by suitably changing the environment, thereby 
renderirig it an important element for the nanobio elec- 
tronics 

An important finding that the overlapping of tt or- 
bitals of the stacked base pairs can create a 7r-pathway 
for charge migration in DNA, has over the years inspired 
many research groups to investigate the electronic prop- 
erties of this important biomolecule 0, |3, H, H, [1,13, H, 
H, [H, [nl, ■ Two independent possible directions can 
be discerned from various experiments, reported as yet: 



Ji3, [S [H, [3 and 

It is now more or 



the charge migration in DNA [4, 
its conductivity [H, HI, [O, [mL 
less established that charge migration in DNA occurs 
thr oug h superexchange tunneling and charge hopping 
d, (Umi. However, the experimental data related to 
DNA conductivity remains unclear: DNA molecules ex- 
hibit a wide range of behaviors, from insulator to metal- 
lic [H, [H, [13, H, HI] ■ Only recently, this issue has re- 
ceived a new twist, i.e., the humidity is recognized as 
an important factor controlling both DNA conductivity 
[H, HE m and DNA magnetic properties [H, [1^. It 
has been observed that the DNA conductivity can in- 
crease exponentially by up to 10^ times with rising hu- 
midity [12, mull. In Ref. [li], participation of the DNA 
molecule in charge transport was verified by the high re- 
sistance of the environment, which exceeds up to 100 



times the resistance of DNA itself. The origin of DNA 
conductivity and its enhancement with humidity was not 
yet clearly understood. The current interpretation of this 
phenomena [1^ rests on the change of DNA permittivity 
and therefore, the DNA conductivity due to adsorption 
of the water molecules on the DNA skeleton. 

For DNA conductivity, the electronic interaction of 
the nearest-neighbor base pairs is the most important is- 
sue, that was extensively investigated by many research 
groups d, [1, [H, [13, [11] . The parameters thus obtained 
were subsequently used for simulations of the charge 
transfer in DNA within different approaches, such as the 
tight-binding Hamiltonian [ij, [H, [13] , the system of ki- 
netic equations [H, [ll| and the polaron model d, [3l| . In- 
terestingly, when the superexchange tunneling and hop- 
ping were taken as two main transfer mechanisms, the ex- 
perimentally observed features for hole migration [1, [l^] 
were indeed reproduced by the theories. However, these 
models fail to explain the diverse behaviors of DNA con- 
ductivity observed in the experiments. 

From our point of view, the inconsistency of theory 
and experiment on DNA conductivity lies with the re- 
ported theoretical approaches, where the charge transfer 
parameters were evaluated via th e q uantum-chemistry 
methods in vacuum [1, [1, [H, [13, [13- Moreover, most 
often this evaluation was performed for a single nucle- 
obase in vacuum [I,[l3,[l8'] that disregards the significant 
shift of the highest occupied molecular orbital (HOMO) 
energies due to the interstrand interaction of the nucle- 
obases participating in the pair formation [1, [13] . The 
discrepancy between the computed parameters and the 
experimental estimate is rather large. For example, the 
potential barrier for the charge transfer from G-C to A-T 
has been estimated theoretically to be ^ 0.7 eV [1, [13|, 
while the experimental value is 0.2 eV [ll[. Application 
of the pure electrostatic model to account for the solva- 
tion effect m, [11] has shown a decrease of this potential 
barrier from 0.7 eV to ~ 0.5 eV, still higher than the 
experimental value. Therefore, it was suspected that 
interaction of DNA with water contributes not only to 
the solvation effect but also by changing the nucleobase 
electronic properties due to their interactions with water. 

Of late several papers have reported quantum-chemical 
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calculations of the electronic properties of DNA sur- 
rounded by solvent molecules [H, [SJ, HI] . We have shown 
earlier [S^ that the nucleobases link to water and inter- 
actions between them can change the symmetry of the 
occupied tt orbitals localized on the nucleobases, that sig- 
nificantly shifts their orbital energies. Therefore, the po- 
tential barrier for hole migration from G-C to A-T pairs 
is decreased from 0.7 eV (dehydrated pairs) to 0.123 (hy- 
drated pairs). Moreover, occurrence of unbound tt elec- 
trons was found to be a reason for an increase of DNA 
conductivity with rising humidity [s^l- The unbound tt 
electrons result in breaking of tt bonds, that occurs be- 
cause of redistribution of the electron density from the tt 
bonds toward the nitrogen atoms due to hydrogen bond- 
ing of nitrogen with water molecules. The increase of 
the charge transfer between two base pairs due to contri- 
bution of such unbound electrons (up to 250 times [s^) 
thus found, is still much less than in the experiments 
[13, m, ■ We suggest that the reason for the under- 
estimation of the humidity effect lies in the fact that the 
DNA geometry transformation from B to A forms with 
dehydration of DNA is not taken into account. 

In our present work we report on our investigation of 
the impact of humidity on the orbital interaction and 
charge transfer between two base pairs if they are stacked 
according to the structural parameters of B-DNA and 
A-DNA. The aspects of the orbital interaction of the un- 
bound electrons are considered also for an explanation 
of the magnetic properties of the B-DNA molecule, ex- 
hibiting the paramagnetic behavior in a magnetic field 

Mm 



II. METHODS 

The molecular dynamics simulations for the DNA 
molecule placed into the "water box" have shown that 
the nucleobases are able to make hydrogen bonding with 
the water molecules [s^. For our investigation of the 
electronic properties of canonical base pairs for dehydra- 
tion and hydration, the optimized structures of the A- 
T and G-C base pairs without water molecules and the 
same structures connected by hydrogen bonds to the wa- 
ter molecules were constructed via the quantum-chemical 
methods. In the first stage, the geometries of the A-T 
and G-C base pairs were optimized in vacuum with the 
Jaguar program [s^]- The Becke3-Lee- Yang-Parr func- 
tionals [3^ and the restricted basis set with polarization 
and difi^use functions 6-31++G** were apphed. Next, 
the crystal structure of the hydrated DNA base pairs 
were obtained by placing the water molecules close to 
the nucleobase atoms having the ability to make hydro- 
gen bonding with the water molecules [3^| . The position 
of water molecules were optimized under the condition of 
frozen base pair geometries to save their planarity. The 
crystal structures of the hydrated (A-T) and (G-C) pairs 
thus obtained are shown in Fig. [TJ 

To build the poly-sequences such as (A-T)2 and (G-C)2 
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FIG. 1: The crystal structures: (a) water-(A-T) and (b) 
water-(G-C) geometries. The water molecules are attached 
by hydrogen bonds to the base pairs at positions given in 
Ref. li^. 



and also the mixed sequences, the optimized geometries 
of the hydrated (A-T) and (G-C) base pairs were stacked 
using the parameters of B-DNA and A-DNA. The natu- 
ral bond orbital (NBO) analysis [3§| for these structures 
was performed, based on the electron density calculated 
with the B3LYP/6-31-I— l-G** functional, which was ear- 
lier found to be a better choice for this purpose [13, [13] ■ 
The application of the HF method for determination of 
the electron density distribution provides similar results 
as the NBO analysis [l^] except for weaker interaction 
of the natural bond orbitals. The eigenfunctions are 
the natural bond orbitals built within the NBO analysis 
based on the input atomic orbital basis set obtained from 
the DFT. The advantage of the NBO analysis is that the 
Pauli exclusion principle is applied not only at the inner 
nodes that preserve orthogonality of the two electrons in 
the same orbital, but also at the outer nodes preserving 
the interatomic orthogonality. Within the NBO analysis 
|,39i] the charge transfer between the nucleobases Q da bas 
been found as the sum of the Vli il* charge transfer 
between the donor orbital Vli belonging to one nucleobase 
to the acceptor orbital il* on another nucleobase as 



(1) 



where qi is the donor orbital occupancy, ei,ej are the 
orbital energies, Ftj is the off-diagonal element. The 
charge occupancy transfer is considered for stabilizing 
the orbital interaction, i.e. when the second order inter- 
action energy AEij = —2Ff^/{ei — ej) is characterized 
by the positive sign. The electronic coupling between 
the nucleobases has been estimated for the natural bond 
orbitals set as 

vda = ^^'^ = E - ^^■) ^ E Q^-i^'^ - ^j-)/?- 



2J 



(2) 

To estimate the tt — tt* charge transfer in the DNA 
molecule, we took into account the charge exchange only 
between the tt orbitals. 



III. TT-ORBITAL INTERACTION 

It is well known that charge transfer in DNA occurs 
due to the overlapping of the tt orbitals of the nearest- 
neighbor nucleobases. The efficiency of charge transfer is 
determined by the coupling of these tt orbitals and the 
difference of their energies. Thus, for two interacting base 
pairs, the transfer of a tt electron between stacked nucle- 
obases is usually attributed to the interaction of HOMO 
and HOMO-1 orbitals (located on purines) and HOMO-2 
and HOMO-3 (located on pyrimidines). The more effi- 
cient charge transfer is expected for two molecular or- 
bitals located on the purines. 

In this section, we analyze the interaction of the 
HOMO and HOMO-1 orbitals in the system of two 
stacked base pairs. The efficiency of such interactions 
directly depends on the location of the base pairs relat- 
ing to each other. It is known that sufficiently wet DNA 
is characterized by a B-form, while dehydration changes 
the geometry of the DNA molecule to the A-form. The 
distinction of A- and B-forms lies mostly in the helical 
twist e, which for A-DNA is 32° and for B-DNA is 36°, 
and in a dislocation of the base pairs from a helix axis A, 
which for A-DNA is 4.5 Aand for B-DNA is -1.0 A[4l|. 
These parameters are important to describe a location of 
a base pair in the system of the two stacked pairs. The 
rest of the parameters, such as base pair tilt 7 and ax- 
ial rise per nucleotide d, which are linearly dependent on 
each other [4l[, do not significantly change this location. 
The geometry of the (G-C)2 sequences, where base pairs 
are stacked with structural parameters of B-DNA and A- 
DNA, are presented in Fig. [2l As the nucleobases in the 
B-form are located on top of each other (the displacement 
from the helix axis A is close to zero) , the most important 
contribution to the base pair interaction comes from the 
intrastrand part. For the A-form, the dislocation of the 
base pair from the helix axis causes both the intra- and 
interstrand charge exchange between the nucleobases. 

Here we study the influence of the base pair location 
in a system of two stacked base pairs, such as (A-T)2 and 
(G-C)2 sequences. The shift of the 7r-orbital energies due 
to the interaction between two stacked base pairs can be 
determined from [i^l 

evi ~ e^J ~ H12 + (e^o - Hi2)Si2 (3) 

e'4>2 ~ (^i>° + H12 - (e^o -I- Hi2)Si2 (4) 

where e^o ^ and e^^ ^ are the 7r-orbital energies respec- 
tively before interaction and after interaction, H12 = 
(-01 I h'^f^ I ?/'2) is the intrinsic interaction integral and 
h*^^^ is the effective core potential. 

In the system of (A-T)2 and (G-C)2 sequences the in- 
teraction of filled TT orbitals is repulsive and therefore, 
the base pairs tend to repeal each other [i^l. This pro- 
cess is dominant when filled orbitals are orthogonalized, 
i.e. when the tilt angle 7=0° and the twist angle 0=0°. 
In this case, because the wavefunctions ij^i and ip2 are 
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FIG. 2: The geometry of the (G-C)2 sequences where the base 
pairs stacked with parameters characterized a) for the B-DNA 
molecule and b) for the A-DNA molecule. 



interacting in such a way that their lobes of the same 
sign overlap (V'l and "02 are orthogonal orbitals, i.e. m- 
phase interaction), this leads to a decrease of the charge 
occupancy transfer Qda between the intrastrand nucle- 
obases and the Pauli repulsion becomes dominant. The 
result of such an interaction is a shift of the 7r-orbital en- 
ergies e^j and e^^ by H12, while the contribution of the 
(e^,o ± Hi2)Si2 term is approximately zero. The shift of 
the HOMO and HOMO-1 orbital energies after perturba- 
tion is presented in Fig. [3] for the (A-T)2 sequence and 
in Fig.|4]for the (G-C)2 sequence. For a twist angle of 0° 
the Tpi wavefunctions are almost equally delocalized over 
two purines and the splitting of the HOMO and HOMO-1 
orbitals is large ~ 2Hi2. 

However, in the A- and B-DNA molecule the pair bases 
are stacked with a helical twist and additionally they are 
shifted from the helix axis. Therefore, the tt electron 
density of the ipi and '02 wavefunctions are shifted such 
that it destroys the wavefunction orthogonality and re- 
sults in a rise of the intrastrand charge occupancy trans- 
fer Qda between the stacked base pairs. Depending on 
the efficiency of the charge occupancy transfer, which is 
determined by the orbital overlap S12 and symmetry of 
the interacting orbitals, the orbital energies e^^ and e^^ 
can be significantly shifted in comparison to the orthog- 
onal case (see the case of 6'=:0° in Figs. [3]- H]). Since the 
TT-electron clouds are shifted, the charge occupancy trans- 
fer is significant even when two orbitals 0i and 02 are 
in-phase (see dehydrated B-DNA in Fig. [5] (a)). There- 
fore, for the B-DNA and the A-DNA, the wavefunctions 
"01 and 02 are strongly localized on one nucleobase and 
sphtting of their HOMO and HOMO-1 orbitals is de- 
creased relating to orthogonal case due to the contribu- 
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b) Hydrated A-T 
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FIG. 3: The shift of the HOMO and HOMO-l energies and 
change of the tpi ^^'^ V'z wavefunct ions' symmetry due to 
the stacking of the two dehydrated A-T pairs (a) (e^o ) or 

1(2) 

two hydrated A-T pairs (b) into (A-T) 2 sequences with dif- 
ferent parameters: 1) two parallel base pairs (7=0°, d=3.4 A) 
stacked with twist angle of 6—0° ; 2) two base pairs staked with 
parameters corresponding to B-DNA (7=0°, d—SA A, 9=36°) 
and 3) with parameters corresponding to A-DNA (7=10°, 
d fa 3.0 A, 61=32°). All energies for the (A-T)2 sequences 
are calculated with respect to HOMO energy of the single 
dehydrated A-T pair (e^o ^ ). 



tion from the (e^o ± Hi2)Si2 term in Eqs. [3}-|4l In our 
earlier works [s^ l , the humidity was found to change the 
electron density distribution and symmetry of the tt or- 
bitals for the single G-C and A-T base pairs. Here, we 
observe that for the (A-T)2 and (G-C)2 sequence in the 
B-form, hydration reverses the interaction of the Tpi and 
-02 wavefunctions. For example, for the (A-T)2 sequence, 
the hydration changes interaction from out-of-phase to 
in-phase, while for the (G-C)2 sequence from in-phase 
to out-of-phase. For the A-DNA the humidity changes 
the symmetry of the '0i and 1^2 wavefunctions as well, 
but preserve the type of their interaction, i.e. the in- 
teraction of the and '02 wavefunctions always remains 
out-of-phase for the (A-T) 2 sequence and in-phase for the 
(G-C)2 sequence. 

The symmetry of the tpi and 1P2 wavefunctions local- 
ized on purines shows an opposite behavior for the (A- 
T)2 and (G-C)2 sequences staked with parameters of the 
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FIG. 4: The shift of the HOMO and HOMO-l energies and 

change of the tpi and tp2 wavefunction's symmetry due to the 

stacking of the two dehydrated G-C pairs (a) (e^o ) or two 

1 (2) 

hydrated G-C pairs (b) into (G-C)2 sequences with differ- 
ent parameters (for parameters see points 1-3 in caption for 
Fig. |3}. All energies for the (G-C) 2 sequences are calculated 
with respect to the HOMO energy of the single dehydrated 
G-C pair (e^o ). 



B- and A-forms. This was found to be a result of the 
interaction of purines with pyrimidine within the DNA 
pairs. For two stacked dehydrated guanines (G)2 and two 
stacked dehydrated adenines (A)2 only out-of-phase in- 
teractions of their ipi and ■02 wavefunctions correspond- 
ing to HOMO and HOMO-l orbitals are observed for 
both A- and B-forms in dehydrated case. Therefore, the 
stacked guanines and the adenines originally are charac- 
terized by the same symmetry of the ipi and ip2 wavefunc- 
tions. Their properties are changed due to the contribu- 
tion of the pyrimidine bases participating in the base pair 
formation. Thus, for the (G-C)2 sequences, the contribu- 
tion of the TT orbital symmetry from cytosine switch the 
interaction of the 0i and -02 wavefunctions to in-phase 
type, while the contribution of the n orbital symmetry 
from thymine for the (A-T)2 sequence do not provide 
such a phase reversal. Similar to the occupied orbitals, 
the (A-T)2 and (G-C)2 sequences have shown the oppo- 
site type of interaction of the lowest unoccupied molecu- 
lar orbitals (LUMOs), such as LUMO and LUMO-1. 

We thus conclude that the transport properties of the 
poly(dG)-poly(dC) and poly(dA)-poly(dT) will be quite 



5 



different due to the difference of their electronic struc- 
tures, namely due to different symmetry of their inter- 
acting HOMOs and LUMOs. This implies that experi- 
mentally observed [lB| different properties of poly(dG)- 
poly(dC) and poly(dA)-poly(dT) molecules showing re- 
spectively the p-type and n-type conductance can be re- 
lated to their intrinsic electronic properties. 



IV. TT CHARGE TRANSFER 

Our next step is the quantitative characterization of 
the charge occupancy transfer between stacked base pairs 
performed with the NBO analysis for DNA of B and 
A forms with different levels of humidity. For the hy- 
drated (A-T)2 and (G-C)2 sequences in the B-form, the 
intrastrand interaction between the stacked pairs and in- 
teraction of the nucleobases with water molecules drasti- 
cally change the electron density distribution and popu- 
lation of the TT orbitals over the base pairs. As a result, 
the covalent structure of the nucleobases is converted to 
a structure with separated charge (ionic), that has been 
obtained by the NBO analysis In particular, for 

adenine and guanine the tt electron density between the 
C4 and C5 atoms is shifted. Therefore, instead of a dou- 
ble bonding to the C5 atom, the C4 atom donates an 
electron to form a double bond with the iV9 atom, which 
contributes one tt electron from the lone pair. Conse- 
quently, the C5 atom has only three bonds and carries a 
negative charge (one tt electron not locked up to the cova- 
lent bond) , while the N9 atom has four bonds (one miss- 
ing TT electron) and carries a positive charge. The charge 
separation occurs for the cytosine in a similar manner, 
where the N3 atom carries a negative charge and the 
NA atom carries a positive charge. Interestingly, for the 
(A-T)2 sequence, the adenines are converted to the ionic 
structure already for the dehydrated case, while thymine 
always has a covalent structure. The origin of this be- 
havior is the low weighting of the covalent structure of 
adenine and the high weighting of thymine [40| . 

The alteration of the double bond pattern due to con- 
version of the covalent structure to the structure with 
separated charges changes the orbital perturbation and 
the intrastrand charge occupancy transfer between the 
nucleobases within each base pair. Therefore, the inter- 
action of guanine and cytosine within the hydrated G-C 
pair (7=0) is significantly increased in comparison to that 
for the dehydrated G-C pair. Their interstrand interac- 
tions are enhanced because of the change of the covalent 
structure of guanine and cytosine to the ionic structure 
and consequently to lowering the energy gap between 
their tt orbitals. In the covalent structure of the de- 
hydrated G-C pair the charge transfer (Qq^c^O. 00043 
e) occurs mostly between the iVl lone pair (guanine) 
and the tt* orbital of the NS-CA bond (cytosine), where 
the energy gap is eAri(G) - eAf3-C4(C) = 5.98 eV. In the 
structure with separated charge, the charge transfer oc- 
curs mostly between the lone pair of the atom (gua- 



nine) and unbound tt electron on the 7V3 atom (cyto- 
sine). The energy gap between them is decreased to 
^Ni{G) ~^N3(C) = 1-9 sV and the amount of charge trans- 
fer is (5g^c=0.0054 e. The electronic coupling between 
guanine and cytosine is estimated to be Vg^c=0.0065 
eV for the hydrated case against Vg^c=0.0016 eV for 
the dehydrated. Unlike in dry DNA in wet DNA the 
interstrand charge transfer within the G-C pair can ac- 
tually make a large contribution to the DNA transport 
properties because of a stronger electronic coupling be- 
tween guanine and cytosine and a lower energy gap be- 
tween them. The enhancement of the interaction of the 
thymine and adenine within a A-T pair due to hydration 
is found to be very small. 

Within the ionic structure of both guanine and ade- 
nine, the unbound tt electrons are localized on the C5 
atom. However, the participation of these tt electrons 
in the intrastrand charge transfer within the (G-C)2 and 
(A-T) 2 sequences is completely different for these purines 
due to the symmetry of the interacting orbitals HOMO 
and HOMO-1 observed in the previous section (in-phase 
and out-of-phase). For the (G-C)2 sequence in B form, 
the "01 and '02 wavefunctions are out-of-phase and there- 
fore, the charge transfer between two guanines is large 
because of the permitted direct charge exchange between 
the C5 atoms (Qc5(Gi)^C5(G2)) with participation of the 
unbound tt electrons. The overlap of these unbound elec- 
trons for the B-DNA is significant (the overlap matrix 
'S'c5(Gi)^ C5(G2)— 0-0245) and energy gap for transfer is 
rather small ec5(Gi) ~ CC5(G2) = 0.33 eV. All of these 
provide a charge occupancy transfer from one guanine to 
another with participation of the unbound electrons to 
be <3c5(Gi)^C5(G2)— 0.145 e, that is a dominant charge 
transfer channel for the (G-C)2 sequence. Because the 
unbound tt electrons on the C5 atoms are strongly lo- 
calized, their orbital overlapping is sensitive to the twist 
angle. Therefore, the Qc'5(Gi)^C5(G2) charge transfer ex- 
ist only in the confined range of twist angle from 0=33° 
to 6=44°. Our computational results for the charge oc- 
cupancy transfer between guanines and cytosines within 
the (G-C)2 sequence are presented in Table HI Dehydra- 
tion of the (G-C)2 sequence in B-form leads to a major 
decrease of the charge occupancy transfer and change 
of its trend as result of application of the symmetry 
rule (see dehydrated B-DNA in Fig. 2]). For the in- 
trastrand cytosines in the (G-C)2 sequences, the domi- 
nance of the 7V3-C4(Ci) C4-7V4(C2) interaction in the 
dehydrated case is switched to the prevalence of the iV3- 
C4(Ci) ^iV3-C4(C2) interaction in the hydrated case, 
that significantly increases the charge transfer between 
the hydrated bases. However, contribution of the cyto- 
sine in the charge occupancy transfer between stacked 
G-C pairs is quit small and can be neglected for both 
B and A-forms. For the (G-C)2 sequence stacked with 
parameters corresponding to the A-DNA, their HOMO 
and HOMO-1 are in-phase that invalidates the charge 
transfer between natural bond orbitals created by un- 
bound electrons and decrease substantially the whole in- 
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TABLE I: The intrastrand charge transfer between guanines 
Gi G2 and cytosines Gi — > G2 within (G-C)2 sequences 
under different conditions. 





Qda, e 


Qad, e 


Q ~ Qda — Qad, e 




Gi ^ 


— » G2 




hydrated B-DNA 


0.1630 


0.0063 


0.1567 


dehydrated B-DNA 


0.0030 


0.0036 


-0.0007 


hydrated A-DNA 


0.0062 


0.0140 


-0.0078 


dehydrated A-DNA 


0.0017 


0.0015 


0.00017 




Ci ^ 


->C2 




hydrated B-DNA 


0.0034 


0.0030 


0.0003 


dehydrated B-DNA 


0.0003 


0.0010 


-0.0007 


hydrated A-DNA 


0.0100 


0.0069 


0.0033 


dehydrated A-DNA 


0.0042 


0.0046 


-0.0004 



trastrand charge transfer between guanines. As was al- 
ready mentioned above, sufficiently wet DNA is charac- 
terized by the B-form and dry DNA is of A-form, and 
the poly(dG)-poly(dC) conductivity with humidity is en- 
hanced 10^-10^ times [22, 23, 24]. The comparison of the 
magnitude Q — Qda — Qad computed for the hydrated 
B-DNA(wet DNA) (see Table IJ) and dehydrated A-DNA 
is calculated to be ^ 10'^ that is in excellent agreement 
with the experimental data [l^]. As was expected, the 
charge occupancy transfer observed between interstrand 
guanine and cytosine belonging to different base pairs is 
quite large {Qq-^_,c2 =0.0322e) for the A-form in compar- 
ison to that for the B-DNA structure ((5gi^C2=0-00136 
e). 

For adenines within the (A-T)2 sequences stacked ac- 
cording to parameters of the B-form, the overlap of 
the unbound tt electrons is S'c5(Ai)-*C5(A2) = —0.0265, 
that is similar to that for the guanines, but ipi 
and '02 wavefunctions are in-phase which restricts the 
C5(Ai) ^C5(A2) charge transfer. However, the un- 
bound TT electrons participate in the intrastrand charge 
exchange with tt* orbitals on the C4-N9{A2) bond, but 
the efficiency is not high, (3c5(Ai)^C4-A'9(A2)=0-0223 e, 
because of a large energy gap ec5(Ai) - ec4-JV9(A2) = 1-9 
eV. Our computational results for the (A-T)2 sequences 
are presented in Table [Tll For adenines within the hy- 
drated (A-T)2 sequences of A-form, the HOMO orbitals 
are out-of-phase and C5(Ai) C5(A2) charge trans- 
fer is permitted. However, the overlap matrix of this 
transfer is small, iS'c5(Ai)^c5(A2) = 0.0130 and the mag- 
nitude of the off-diagonal element -Fc5(Ai)->C5(A2) is 10 
times smaller than that for B-DNA. Therefore, the con- 
tribution of the charge exchange between unbound elec- 
trons in (A-T)2 sequences of the A-form is less. Because 
thymine structure remains covalent independent of hy- 
dration, the charge transfer between thymines for the A 
and B-form is slightly modified by hydration only due to 
a change of the (e^ — Sj) energy gap. Finally, an increase 
of the TT — TT* charge transfer for the (A-T)2 sequences 
is found to be only ~ 40 times due to the hydration, 
while the experimental value was 10^ times This 
discrepancy is related to the applicability of the simu- 



TABLE II: The intrastrand charge transfer between adenines 
Ai A2 and thymines Ti — > T2 within (A-T)2 sequences 
under different conditions. 





Qda, 


Qad, e 


Q = Qda — Qad, e 




Ai - 


A2 




hydrated B-DNA 


0.0269 


0.0052 


0.0218 


dehydrated B-DNA 


0.0040 


0.0026 


0.0014 


hydrated A-DNA 


0.0053 


0.0140 


-0.0087 


dehydrated A-DNA 


0.0025 


0.0055 


-0.0030 


Ti ^T2 


hydrated B-DNA 


0.0020 


0.0010 


0.0010 


dehydrated B-DNA 


0.0018 


0.0011 


0.0007 


hydrated A-DNA 


0.0056 


0.0033 


0.0023 


dehydrated A-DNA 


0.0057 


0.0033 


0.0024 



lated results only for low temperature range, where or- 
bital symmetry of the unbound electrons is preserved, i.e. 
when the energy required to change the orbital symme- 
try is AEs <C ksT/q. We estimated the magnitude of 
AEs for adenines within the (A-T)2 sequence as the dif- 
ference between the ec5(Ai) — ec5(A2) for unbound elec- 
trons when the wavefunctions are out-of-phase and in- 
phase, which was found to be Ai?5=0.12 eV. There- 
fore, already at a low temperature (T <100 K) unbound 
TT interacting electrons have enough energy to switch 
the wavefunctions from out-of-phase to in-phase, which 
opens up the opportunity for a direct C5(Ai) — !-C5(A2) 
charge transfer in the (A-T)2 sequences of the B-form. 
As a result, at room temperature the charge transfer 
with participation of the unbound electrons in B-DNA 
is possible not only for the (G-C)2 sequence but also 
for the (A-T)2 sequence. Therefore, at room tempera- 
ture the hydration of DNA should provide an increase of 
the DNA conductance approximately as the same magni- 
tude for the poly(dG)-poly(dC) and poly(dA)-poly(dT) 
chains (up to 10'' times), because the orbital overlap 
of nearest-neighbor unbound electrons for stacked gua- 
nines and stacked adenines is the same {Sc5(Gi{Ai))^ 
C5(G2(A2)) ~ ±0.025) and orbital symmetry is not pre- 
served at room temperature. 

In summary, we found that the occurrence of the un- 
bound TT electrons activated by humidity suppresses the 
DNA band gap from 8.0 eV for the dehydrated DNA 
to ~ 3.0 eV for the hydrated DNA. This value of the 
band gap was often observed in the experiments [H, Q . 
Therefore, the occurrence of the unbound tt electrons 
suppressing the band gap of wet DNA is the main factor 
for high B-DNA conductivity. 

V. DNA MAGNETISM 

A noninvasive means to investigate the intrinsic elec- 
tronic properties of DNA was recently applied to wet 
and dry A-DNA by measuring the magnetization of this 
molecule at different humidity levels [25|. The paramag- 
netic behavior of A-DNA in the B-form was observed at 
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low temperature, while dry A-form was shown to be dia- 
magnetic for a wide temperature range. The solvent of 
distilled H2O itself is paramagnetic due to the contribu- 
tion from the O2 molecules, but the measured magnitude 
of magnetization of distilled water was not large enough 
in comparison to that for wet A-DNA [26|, which proves 
the origin of the DNA paramagnetism from the intrinsic 
DNA properties. The orbital motions of free electrons 
have been attributed to be the reason of this param- 
agnetic effect, while the spin contribution to paramag- 
netism was claimed to be unlikely. However, the interac- 
tion of nearest-neighbor unbound tt electrons created by 
interaction of the base pairs with water molecules can be 
a source of spin paramagnetism, because of possible pair- 
ing of the unbound electrons with parallel or anti-parallel 
spins. 



We have shown in the previous section that each DNA 
base pairs linked to water molecules has unbound tt elec- 
trons, whose interaction is not strong because of the large 
separation (3.4 A). The consideration of the Pauli exclu- 
sion principle for outer nodes (interatomic orthogonality) 
in the NBO analysis allows us to analyze the spin state of 
unbound electrons. Originally, the phase of the interact- 
ing wavefunctions is taken from the standard calculation 
B3LYP/6-31-I-I-G** within the restricted basis set. If 
the wavefunctions of the HOMO and HOMO-1 orbitals 
obtained within the DFT method are in-phase, then the 
unbound electrons obtained within the NBO analysis are 
also in-phase (see Figures [3] and |4]), i.e., paired with 
parallel spin. Their natural orbitals are closer in energy 
because of absence of exchange repulsion term K in the 
interaction energy of two electrons. For the two unbound 
electrons paired with anti-parallel spins (out-of-phase in- 
teraction of HOMO wavefunctions) their orbital energy 
splitting is quite large. The energy diagram of the un- 
bound electrons in the guanines within the (G-C)2 se- 
quence, in adenines within the (A-T)2 sequence and in 
cytosines within the (G-C)2 sequence are presented in 
Fig. [5] a), b) and c) respectively for hydrated sequences 
in B and A-forms. Two stacked cytosines are found to be 
source of paramagnetic behavior because their paired un- 
bound electrons have parallel spin, while for the stacked 
adenines and guanines the interaction of the unbound 
electrons largely depend on the DNA geometry. However, 
the stacked adenines as are also expected to contribute 
to the B-DNA paramagnetism. 



Therefore, with regards to the diamagnetic and para- 
magnetic properties of DNA, dry DNA is expected to be 
always diamagnetic because of the absence of unbound 
electrons, while wet DNA should be paramagnetic and 
the efficiency of paramagnetism is directly dependent on 
the DNA sequences. The control of the spin of the un- 
bound electrons by the magnetic field would open up real 
opportunities for using DNA in spintronics [45|. 



a) Guanine, K=0.15 eV 

Ec5. eV g.DNA 

4 



A-DNA 



-3.25 - 



-3.50 - 



-3.75 H 



t 



3.34 eV 
0.35 eV 

3.69 eV 



-3.32 eV 
- 0.2 eV 
-3.52 eV 



b) Adenine, K=0.12 eV 

Ec5. eV 

-3.00 H B-DNA 



-3.25 - 



-3.50 H 



0.07 eV 

-3.32 eV 
-3.39 eV 



A-DNA 
4 



-3.13 eV 
- 0.19 eV 
-3.32 eV 



c) Cytosine 

En3, eV B_[5N^ 



A-DNA 



-5.40 - 



-5.50 



-5.60 H 



- 4- 



-5.44 eV 
-5.49 eV 
0.05 eV 



0.04 eV 
-5.55 eV 
-5.59 eV 



FIG. 5: The energetics of the natural bond orbitals Q.i corre- 
sponding to the valence unbound electrons formed within the 
nucleobases due to their interaction with water a) for gua- 
nines within (G-C)2 sequence, b) for adenines within (A-T)2 
sequence and c) for cytosines within (G-C)2 sequence. K is 
the exchange repulsion term calculated as difference of energy 
AiJ for paired electrons with anti-parallel spins and parallel 
spins. 



VI. CONCLUSION 

Our findings that hydration of DNA activates the oc- 
currence of the unbound tt electrons, which can easily 
participate in conductance because of narrowing of the 
band gap to ~ 3.0 eV in comparison to ^ 8.0 eV for dehy- 
drated DNA, allow us to characterize the DNA molecule 
as a narrow band gap semiconductor. The main factors 
enhancing the appearance of these unbound electrons are 
linking of the nucleobases with water molecules and in- 
trastrand interactions between neighboring nucleobases. 
The main contribution to the conductivity comes from 
the intrastrand transfer of unbound tt electrons between 
purines because of quite strong overlapping of their or- 
bitals (S'c5(Gi(Ai))^C5(G2(A2)) ~ ± 0.025) and a smaU 
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energy gap ec5(Gi(Ai)) - ec5(G2(A2)) = 0.33 eV. There- 
fore, at room temperature sufficient hydration of DNA 
can lead to an increase of the conductivity by up to 10'^ 
times, but at low temperature the conductance proper- 
ties of the poly(dG)-poly(dC) and poly(dA)-poly(dT) se- 
quences should be significantly different because of the 
orbital symmetry rule. As an example, because of the 
out-of-phase interaction of the wavefunction of the cor- 
responding TT unbound electrons, the conductivity of the 
poly(dG)-poly(dC) chain is expected to be much larger 
than that for the poly(dA)-poly(dT) structure, where 
these wavefunctions are in-phase. The pairing of the un- 
bound TT electrons from the nearest-neighbor intrastrand 



bases within the hydrated DNA chain can be with anti- 
parallel or parallel spins, that is a source of correspon- 
dent diamagnetic and paramagnetic behavior of DNA in 
a magnetic field. 
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